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It has recently been argued that bubble nucleation in ekpyrotic and cyclic cosmological scenarios 
can lead to unacceptable inhomogeneities unless certain constraints are satisfied. In this paper we 
show that this is not the case. We find that bubble nucleation is completely negligible in realistic 
models. 


The cyclic model of the Universe 1,2 is a radical al¬ 
ternative to standard big bang and inflationary cosmol¬ 
ogy. The cyclic picture proposes that the Universe under¬ 
goes an endless sequence of cosmic epochs which begin 
with expansion from a big bang and end in collapse to 
a big crunch. The connection between the big crunch 
and the ensuing big bang is presumed to occur accord¬ 
ing to the prescription recently proposed with Khoury, 
Ovrut, Seiberg. 3 The expansion part of each cycle has a 
period of radiation and matter domination followed by 
an extended period of cosmic acceleration at low ener¬ 
gies. The long period of acceleration is crucial in estab¬ 
lishing the flat and vacuous initial conditions required 
prior to beginning the contraction phase. The density of 
entropy, black holes and other debris from the previous 
cycle is reduced to nearly zero. Subsequently, the Uni¬ 
verse ceases to accelerate, fluctuations are generated, and 
the Universe heads towards a big crunch in which mat¬ 
ter and energy are regenerated. The key ingredient is the 
restoration of the Universe to a pristine flat vacuum state 
before each big crunch, which ensures that the cycle can 
repeat and that the cyclic solution is an attractor. 

The cyclic model can be described in terms of a scalar 
field <f> moving back and forth in an effective potential 
V{<t>) (see Fig. 1 ). Although the model can be de¬ 
scribed in terms of an ordinary quantum held in four 
dimensions, the model is strongly motivated by ideas de¬ 
veloped as part of the ekpyrotic scenario, 4 drawing on 
the braneworld picture and string theory. The scalar 
held represents the modulus held that determines the 
distance between branes. Moving from some positive 
value to minus infinity and back along the potential cor¬ 
responds to a pair of orbifold branes being drawn to¬ 
gether by an interbrane force, colliding (corresponding 
to cf) — > — oo), bouncing and rebounding. The transi¬ 
tion between big crunch and big bang corresponds to the 
bounce. 3 The accelerating phase begins when the branes 
are maximally separated (cf> = <f> o) once the matter den¬ 
sity falls below the positive interbrane potential energy 
density. As the branes are drawn together, the accel¬ 
eration rate decreases until, at <f> = <f>i, the expansion 
begins to decelerate. As the branes draw yet further 
together, the potential becomes negative (<j> < 0), the 
expansion stops (</> = <(>2) and the universe begins to con¬ 


tract. Fluctuations are generated during the contraction 
phase ((/) = <f> 3), before the potential reaches a minimum 
at (p = min■ 


m 



FIG. 1. Plot of the actual analytic potential (solid line) 
described in Refs. 1 and 2, shown to proper scale. For the 
purposes of analyzing bubble nucleation, differences between 
the schematic plots shown in Refs. 1 and 2 and the prop¬ 
erly scaled plot are significant. The axes are expressed in 
Planck units. The value of pq is about + 10 -120 correspond¬ 
ing to the current dark energy (quintessence) energy density. 
The markings (fn refer to points described in the text. The 
dotted line is the piecewise linear approximation assumed by 
Heyl and Loeb. Our central point is that the piece-wise lin¬ 
ear function is an exponentially poor approximation to the 
slope for <j> < (j>\, leading to an exponential overestimate of 
the nucleation rate. If the kink is moved to cj>i, the piece- 
wise approximation is reasonable, but the decay rate is then 
negligible. 

Recently, Heyl and Loeb 5 have argued that bubble nu¬ 
cleation during the accelerating and collapse phase, near 
<j> = <f> 1, for example, can spoil the homogeneity of the 
universe unless the potential obeys certain constraints. 
They have focused on tunneling about </> 1, where the the 
kinetic and potential energy become comparable as the 
field rolls to the left and where the accelerating expansion 
phase ends. Bubbles form whose interior corresponds to 
values of <f> far down the potential. Once formed, the 
bubbles grow at the speed of light.) These bubbles de¬ 
stroy the homogeneity essential to the scenario. Their 
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point is that, although the bubbles are rare events, there 
remain more than a billion years from when passes (f>\ 
until it reaches the bottom of the potential well. They 
estimate that this is sufficient time to have a significant 
amount of bubble nucleation, unless the potential well is 
not very deep. 

In this paper we show that the approximations used by 
Heyl and Loeb, although motivated by schematic figures 
of the potentials in Refs. 1 and 2, are not valid for the 
actual analytic forms, as presented for either ekpyrotic 
or cyclic models. In particular, they estimate the bubble 
nucleation rate assuming a potential which can be ap¬ 
proximated as two linear segments, one with slope zero 
and one with non-zero slope (Fig. 1). However, the slope 
of the potentials used in the ekpyrotic and cyclic models 
is lO 100 or so times smaller than in their piece-wise lin¬ 
ear approximation for cf >4 <C (f> < 4>\ . This error is further 
amplified in the computation of the bubble nucleation 
rate, which is exponentially sensitive to the slope. The 
smaller is the slope, the more suppressed is the nucle¬ 
ation rate. For the potentials actually used, we will show 
that the bubble nucleation rate is completely negligible 
compared to the classical rolling rate. 

The explicit model for V (</>) which was discussed in the 
context of the cyclic model is of the form : 1,2 

V(<j>)=Vo(l-e- c *)F(ci>), (1) 

where, without loss of generality, we have shifted </> so 
that the zero of the potential occurs at (f> = 0. The 
function F((f>) is introduced to represent the vanishing of 
non-perturbative effects described above: F(cj>) turns off 
the potential rapidly as </> goes below <fimin, but it ap¬ 
proaches one for </> > 4> m in- For example, F{(f >) might 
be proportional to e~ l ^ 9a or e -1 / 9s (for illustration, we 
will assume the latter form), where g s oc e 1 ^ for 7 > 0 . 
We will choose c = 10 and 7 « 1/8 so that the poten¬ 
tial minimum corresponds to V((j> m i n ) = M 0 4 = 0(1 ) in 
Planck units, a choice which Heyl and Loeb claim is ruled 
out by the nucleation constraint. The constant Vq is set 
roughly equal to the vacuum energy observed in today’s 
Universe, of order 10 “ 120 in Planck units. For large c, 
this potential has V"/V 7C 1 for <f> 1 and V"/V >> 1 

for (fimin < <f> < 0. These two regions account for cos¬ 
mic acceleration and for ekpyrotic production of density 
perturbations, respectively . 4,6 More precisely, the end of 
acceleration is at </> 1 « 0.35. The rest of the markers in 
Figure 1 are scaled to the same parameters. The cur¬ 
rent value (circle) corresponds to tens of trillions of years 
before the big crunch. In the ekpyrotic region, the con¬ 
stant term is irrelevant and V may be approximated by 
—Vo e~ c The scaling solution over this range is 

a(t)=\t\r, V = -V 0 e- C * =- P{1 ~ 3P \ p=|. 

( 2 ) 


Here we chosen t = 0 to be the bounce. These expressions 
are used to obtain the time estimates discussed in the 
paper (and this agrees with Heyl and Loeb). 

Heyl and Loeb focus on (j) near <j> 1 , where the scaling 
solution first becomes valid. The Universe is expand¬ 
ing but at a decelerating rate. As the field roll towards 
more negative values, the potential energy becomes in¬ 
creasingly negative, the expansion stops and contraction 
begins. They approximate the potential as 

V(4>) = -M 0 4 Q(</>1 - <t>) (3) 

tymin 


where ©(</>) is the Heaviside step function. This amounts 
to approximating V(<j)) =0 for </>></1 and U(</>) = 
K ■ (<j> — </>! ) for (f> < <f>i , where K is a constant chosen to 
be the slope of the line that connects </> 1 to <t> rn in • We have 
superposed this approximate form in Figure 1. The piece- 
wise linear potential was introduced by Lee and Weinberg 
in their analysis of bubble nucleation for a rolling scalar 
field without a barrier . 7 

There are two important aspects to the bubble nucle¬ 
ation problem in the piece-wise linear potential. First, as 
derived by Heyl and Loeb , 5 gravitational effects suppress 
the tunneling altogether unless <f> is within a Planck dis¬ 
tance of <f>i, 77 = (j> — </>i < 0(1). Hence, for values of q J 
far to the right of <j> 1 on the plateau, tunneling past 1 is 
insignificant even though the classical rolling rate is slow. 
Second, as derived by Lee and Weinberg, the bubble in¬ 
terior corresponds to a jump in <j> in which <j> jumps past 
</>i by an amount 77 . So, the closer is <f> to (j> 1 , the less is 
the jump. Combining the two aspects, it is clear that a 
good piece-wise linear approximant to V((/>) to be used to 
estimate tunneling near <j> 1 should match the slope near 
<j> 1 within one Planck distance, A <f> < 1. The behavior of 
V((f>) for values farther away in either direction plays no 
role in the nucleation rate about (pi. 

As is evident from Fig. 1 , the slope of the piecewise 
linear approximation disagrees by an exponential factor 
from the actual slope of the potential for </><</> 1 . Within 
a Planck distance of <j>i , the difference in the actual slope 
and the piecewise approximant is a factor of nearly 10 lo °. 
The difference totally changes the conclusion. The bub¬ 
ble nucleation rate is proportional to 7 


exp 


327T 2 (4> — <pi) 3 
~ 3 K 


(4) 


That is, the exponentially smaller slope of the actual po¬ 
tential is further exponentiated in computing the nucle¬ 
ation rate. Consequently, the nucleation rate is exponen¬ 
tially smaller than the estimate obtained by Heyl and 
Loeb and is completely negligible, as we suggested in our 
paper . 2 

The argument is not complete. We have argued that 
the piecewise linear approximation is poor near <f> = (f> 1 . 
However, one might worry that it is a reasonable ap¬ 
proximation for (j> near ((> 4 , where the potential becomes 
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steeper. Does nucleation become significant as 4> ap¬ 
proaches </> 4 ? First, note that Heyl and Loeb’s result can 
be generalized to show that nucleation is insignificant if 

(<j> - <^4) 3 (<^4 - (pmin ) > Mq (5) 

where — Mq corresponds to the minimum of the potential. 
Here we have chosen Mo in Planck units. In this case, 
where (f >4 — 4> m in ~ 3, we can ignore nucleation unless <f> 
is within a few Planck distances of <f> 4 . The stronger con¬ 
dition, cited above, is that gravitational effects suppress 
nucleation unless <f> is less than one Planck unit from <j> 4 . 
Then, the nucleation rate is comparable to what Heyl 
and Loeb found. However, the time left before the field 
rolls to potential minimum is exponentially smaller. At 
(j> m <f> 1 , the remaining time is about Hq 1 = 10 60 in 
Planck units. For <fr within less than a Planck distance 
from <f> 4 , the remaining time is less than one Planck time. 
There is not enough time left to have significant chance 
to nucleate because here the field is rolling too quickly. 

If one considers tunneling near points between (j>\ and 
(f >0 (today’s value), the tunneling rate is more suppressed 
than at <f>\ because the slope of the potential K is less. 
The rolling time increases as (U ") _1 « 1/cK but the 
tunneling time increases by ~ exp(l/AT), an exponen¬ 
tially greater factor. 

We confess to bearing some responsibility for Heyl and 
Loeb’s mis-application of the piecewise linear approxi¬ 
mation to the potential. We provided in our papers 1,2 a 
schematic drawing of the potential which suggests that 
the slope near <j) = <f> 1 (where there is a billion years be¬ 
fore reaching the minimum) is steeper than is actually 
the case. This is apparent if one plots the actual poten¬ 
tial function described in the text, as has been done in 
Figure 1. 

The bottom line is that bubble nucleation is negligible 
throughout the cyclic solution and provides no signifi¬ 
cant contraint for the potentials we consider in ekpyrotic 
and cyclic models. However, it is good to be aware of 
Heyl and Loeb’s analysis if one considers more general 
potentials. 
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